Self-similarity in a system with a short-time delayed feedback 
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Using the Poincare section technique, we study in detail the dynamical behaviors of delay dif- 
ferential system and find a new type of solutions Si in short-time delay feedback. Our numerical 
results remind us to deny the opinion that there are no complex phenomena in short-time delay 
case. Many similarities between foundamental solution and the new type of solutions are found. We 
demonstrate that the scales of Si increase with exponential growth via i in the direction of fj,, while 
decrease with exponential decays in the direction of x or delay time tn. 
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Optical feedback systems governed by delay differential equations (DDEs) have attracted much attention from both 
■ the applied and the fundamental points of view [1-16]. Generally, the delay-differential system related to optical 
0\| bistable or hybrid optical bistable device is described by 

t-H T , x(t) = -x(t) + f(x(t-t R ), f i), (1) 

> ; 

where x{t) is the dimensionless output of the system at time t, t R is the time delay of the feedback loop, r' is the 
response time of the nonlinear medium, the parameter a is proportional to the intensity of the incident light. In Eq. 
([j]), f(x,fx) is a nonlinear function of x, characterizing the system, e.g. f{x,fi) = /U7r[l — Ccos(a; — xb)] for Ikeda 
f - ^ ■ model Q], f(x, fx) — n[A — (ism 2 (x — Xq)] for Vallee model and f(x, fj,) = fism 2 (x — Xq) for the sine-square model 

o>: i- ^ 

ON ■ The understanding of the Eq. ([1]) up to now can be summarized as follows. The first experimental observation of 
period-doubling bifurcations and chaos in a hybrid bistable device was made by Gibbs et al. |jj following a prediction 
by Ikeda et al. ||. The solution of system, which appears after Hopf bifurcation, evolves through a period doubling 
Tp — » 2Tp ■ ■ ■ — > 2 N Tp, as one parameter is varied. These solutions are called 2 N periodic and the cascade accumulates 
at the Feigenbaum point. These solutions are named fundamental solutions by Ikeda et al. we do so in this paper. 
Later the two groups found that higher- harmonic oscillation states appear successively in the course of transition to 
developed chaos in long-time delayed case (i.e. delay time is longer enough than the response time). These solutions 
coexist and each follows a period-doubling cascade. The oscillation period of these harmonic states are given by 
Tp In, where n is odd integer and Tp is the period of the fundamental solution. As the study of the dimension of the 
chaotic attractor, the behavior of the DDE exhibits high-dimensional chaos [[|. Ikeda and Matsumoto have given an 
estimate of the Lyapunov dimension of attractor for the Ikeda model, and it ranges approximately from 2 to 13 when 
some bifurcation parameter is varied. Recently, some researchers demonstrated that the behavior of quasiperiodicity 
followed the hierarchy of the Farey tree |^-[l0|] and the chaotic itinerancy phenomenon switches among some different 
unstable local chaotic orbits pT| - b| . We reported two new types of solutions found in moderate-time and short-time 
delay regimes [ ps| , which are different from the fundamental solution and the odd harmonic solution. In this paper, 
we study in detail the dynamical behaviors of the new type of solution found in short-time delay regimes. 

This paper is organized as follows: In Sec. II, the numerical methods used in this paper are introduced. By using 
Poincare section technique, we can easily observe the course of bifurcation of DDE, and easily distinguish the new 
type of solutions and the fundamental solution. In Sec. Ill, we demonstrate our numerical results. In the short-time 
delay case, there is a new type of solutions Si(i = 1,2, 3...), which has many similarities comparing with fundamental 
solution. Moreover, these new solutions are alike each other. In Sec. IV we summarize our results and conclude. 
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II. NUMERICAL METHODS 



Measuring the delay time in units of Ir, one can rewrite Eq. (|l|) as 

Tx(t) = -x{t)+f(x{t-l),n), (2) 

where r = r'/tR characterizes the effect of the time delay when r' is fixed. In this paper, we study Eq.(Q) with the 
special feedback function 

f(x,fj,) = 1- fix 2 . (3) 

This feedback function can be considered as the first nonlinear term of the Taylor expansion of the general nonlinear 
function f {x , /i) in the vicinity of a steady state. It should keep the general nonlinear properties of DDE, as shown in 
Refs. [0]l5). 

The Eq. (Q) can be solved numerically and a fourth-order Adam's interpolation is suitable for that. In order to 
trace the evolution of a DDE, one might investigate the evolution curve of the variable x(t) vs the time t. However, 
it is difficult to distinguish different solutions if one only observes the x(t) — t relation. Some of us (Zhao et al.) have 
offered a method in Ref. jio) to represent the solutions of a one- variable DDE by using the Poincare section technique. 
This method has been proved to be a powerful tool in exploring the evolution of bifurcation of DDE. Let us review 
this method briefly. Let x t {9) = x(t + 9), — 1 < 9 < 0, then x t2 (9) is determined by x tl (9) uniquely according to Eq. 
(||), where t\ < £2- We approach the section mapping as follows: choose an appropriate constant x c € R; integrate 
Eq. (J2J) numerically till x(t) > x c and x(t + h) < x c , where h is the length of the integrating step; then proceed a 
simulation procedure to get ti as well as x ti {9) such that x ti (0) = x c . To be simple, we denote x ti (9) as Xi(9) in 
the following discussions. In this way we convert the flow of Eq. (Q) into a mapping which maps the curve Xi{9) 
onto the curve Xi+\(9). We regard this curve-to-curve mapping as the Poincare map of a DDE. A periodic solution 
of Eq. (||) with period T, x(t) = x(t + T), corresponds to a periodic solution of the Poincare map with period N, 
Xi{9) = Xi+N{9), where N is an integer. For practical applications, we can take n discrete points Xi(9j) on the curve 
Xi{9) to represent the solution, where 6j G (—1,0) and j = 1,2, ...,n. Then the curve-to-curve mapping appears as 
a point-to-point mapping in R n . In order to exhibit the bifurcation process, here we usually need a one-dimensional 
mapping representation Xi{9\) with the bifurcation parameter. 

III. RESULTS 

As usually considered, there is no complex phenomenon in short- time delay region since Eq. (^|) will approach a 
normal one-dimensional ordinary differential equation. Our results remind us that this is not the truth. 

A. The bifurcation of fundamental solution 

Before discussing the new type of solutions, let us first review the bifurcation process of the fundamental solution. 
In the long-time delay case (i.e. r is very small), Ikeda et al. Q have shown theoretically that instabilities and chaotic 
behaviors can occur in the system. As t is fixed and /1 is increased, a square-wave solution appears after the Hopf 
bifurcation of a steady state. With further increase of /i, this square wave solution undergoes a square of bifurcation 
with its period doubling itself successively and then becomes chaotic. We define this solution as the fundamental 
solution of the system and marked it as 5*0- 

When t increases (i.e. delay time decreases), the fundamental solution exhibits mirror-similar bifurcation behavior 
as shown in Fig. 1(a). With the continuous increase of t, the period-doubling bifurcation with less and less order 
takes place in the course of bifurcation. At t = 1.13, So undergoes only period-two bifurcation via /1. Fig. 1(b) shows 
a bifurcation diagram just below this value. With further increase r, we can no longer observe the period-doubling 
bifurcation of So, see Fig. l(c)-(f). We regard the regime of r > 1.13 as short-time delay case. In this regime, the 
delay time tn is smaller than the response time r' because r = r'/t/j, and the fundamental solution exhibits only 
one-period limit cycle state. 

B. New type of solutions 

In the short-time delay case, fundamental solution shows no bifurcation and chaos. This is not to say that there is 
no chaos in this system with the varying of \i at a fixed parameter r. In fact, there still exists another chaotic attractor, 
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which locates behind So in the direction of fx, as shown in Fig. l(c)-(f). We marked them as Si(i = 1, 2, 3, ...). In the 
direction of /Lt,every Sj evolves a period-doubling bifurcation. Figure 2 exhibits the evolution courses of fundamental 
solution and Si, which is located on period-one limit cycle state of themselves, respectively. 

From Fig. 1 and Fig. 2, we can easily find that there are not only similarities but also difference among So and Si. 
Firstly, with the increase of r, each Si appears the same bifurcation process as Sq does. Comparing Fig. 1(c) with 
Fig. 1(a), we can find the diagram of Si at r = 3.0, displays the same shape as that of So at r = 0.80. In fact, at 
certain parameter, S2, S3 and S4 et al. also show the similar shape. Secondly, Si has more and more oscillation with 
the increase of the subscript i. Simultaneously, i increases with the increase of r. So has only one peak within one 
period. In contract to So, Si has not only the peak but also another small peak; S2, S3 and S4 have more and more 
peaks within one period respectively, see Fig. 2. Thirdly, with the increase of t, Si appears one by one and Si follows 
So, S2 follows Si in the direction of /x. 

C. The scales of Si via i 

From Fig. 1 one can find that S,; appears continually with the increase of r. In order to find the law which exhibits 
the appearance order of Sj and the scales of Si, we should choose a standard to compare Sj with each other. In this 
paper, we choose the critical values of Tj as the standard, Tj is the value at which the second period-doubling bifurcation 
of the period-1 solution of Si takes place with the decrease of r. In fact one could also choose another standard. For 
this specific choice, the bifurcation diagram of Si appears as the patterns in Fig. 3(a)-(e) respectively. Our numerical 
solutions show that n increases with exponential growth against the increase of i. Figure 4(a) demonstrates the result, 
where the scatters are the values of Ti and the dotted line is the fitting curve which is exponential growth function 
with the increase of i as follows: Ti = Aexp(i/B), where A and B are fitting coefficients.. 

From Fig. 1 one can find the scale of the bifurcation diagrams of Sj along \i increases, while along x decreases 
with the increase of i. Using the standard chosen above, we can measure the length Sfj, of the period-2 solution in the 
direction of /1 at Ti and use it to characterize the scale of Sk along \x. Figure 4(c) shows that Sfi also increases with 
exponential growth via i, and the exponential function is Sfi(i) — Aexp(i/B). On the other hand, by measuring the 
maximum height 5x of the period-2 solution at Ti, we find that Sx exhibits the exponential decay with the increase of 
i: 8x(i) — Aexp(— i/B),see Fig. 4(d). 

IV. CONCLUSION 

In this paper we have studied in detail the dynamics of short-time delay differential system. Our numerical result 
remind us that it is not truth to consider that there is no complex behavior in short-time delay feedback. By using the 
Poincare section technique in DDE ]ic| , a new type of solutions Sj was found in this case and it has many similarities 
as same as fundamental solution in the bifurcation diagrams. We found the law of Si and showed the scales of Sj 
with the increase of i. The scales of Sj increases along fi and t, while decreases along x and delay time . 
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FIG. 1. The bifurcation diagrams of So and Si with the increase of r. 



FIG. 2. The evolution of So and S t versus time, (a) S ,r = 0.3, fj, = 1.50; (b) S u t = 1.5, fx = 8.50; (c) S 2 ,r = 1.5,// = 9.42; 
(d) S 3 , r = 3.5, /i = 31.45; (b) St, r = 40, fi = 2585. 



FIG. 3. Bifurcation diagrams of So and Si at the critical values of n. 



FIG. 4. (a) n, (b) tn, (c) Sfi, and (d) 5x of Si versus i. Here n is as same as figure 3 and tn = r'/r. Sfx is the length 
and Sx is the maximum height of period-2 solution in Fig. 3 respectively. The square scatter is the numerical result and the 
dotted line is the fitting curves, (a) n = 2.026 exp(i/0.817); (b) t R (i) = 1.135 exp(-i/0. 747); (c) Sfi(i) = 1.072 exp(i/0. 616); 
(d) 5x(i) = 0.125 exp(-i/0.458). 
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Fig. 1 
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